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Abstract 

I propose a stochastic SIS and SIRS system to include a Poisson measure term to model anomalies in the 
dynamics. In particular, the positive integrand in the Poisson term is intended to model quarantine. Conditions 
are given for the stability of the disease free equilibrium for both systems. 
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1 Introduction 

Many authors have considered continuous time stochastic epidemiological models [TTJ [TSJ [T7] , however, none of these 
models taken into account anomalies that affect the system. Examples of these events are super-carriers, the amount 
of contacts of individuals in the stadium and in the bars during a sports event, quarantine, etc. A particular example 
is the three waves of influenza pandemic in 1918 [5]. 

For a fixed population of size N, we consider the following SIS model: 

S'(t) = -PS(t)I(t) - »S(t) + ^ + XI(t), 

l'(t) = pS(t)I{t)-{\ + »)I{t), u 



and the SIRS model: 

S'(t) = -0S(t)I(t)+5R(t), 

I'(t) = (3S(t)I(t) ~ XI(t), (2) 
R'{t) = XI(t) - 5R(t). 

Here, S(t), I(t), and R(t), denote the frequencies of the susceptible, infected, and removed respectively, where 
S(t) +I(t) + R(t) = 1. The constant fj, represents the birth and death rate (newborns are assumed to be susceptible), 
A is the recovery rate for the individuals that are infected, /3 represents the average number of contacts per day, and 
5 is the rate for which recovered individuals become susceptible. 

For the SIS model, if j3 > fi + X, then an epidemic will occur, and if /3 < fj, + A then the process will converge to 
the disease-free equilibrium, which means that the disease has disappeared. In other words, if /3 < /i + A then point 

(1, 0) is globally asymptotically stable, and if /3 > fi + A, the point — ? 1 — — ) (the endemic equilibrium) is 

globally asymptotically stable. 

Considering the SIRS model, when /3 < A the disease free equilibrium is globally asymptotically stable. If /? > A 

fX X 1 — — 1 — — \ 
then the endemic equilibrium ^— , 1 — — — j, jj is globally asymptotically stable. 

We now define and give some intuition to the Poisson measure. Take X to be an K-valued Levy process on 
(f2, J 7 , P) with the filtration {•7 r t}teR + , where Tq contains all of the null sets of T. For B 6 £>^K rf \{0}^ and a fixed 

w e O, define AX u (s) := X u (s) - X u (s-) and N u (t, B) := #jfj < s < t : AX u (s) G Bj. This random counting 
measure is known as the Poisson measure since, fixing B, the map t —> A^(£, B) is a Poisson process with intensity 
tv(B) := E N u (t, B) . We call u(-) the intensity measure and it is well known that it is a Levy measure, hence, 



1 



it is a Borel Measure with / I \x\ 2 A l)u(dx) < oo. The measure N u (dt, dx) :— N U! (dt,dx) — dtv{dx) is called 

Jk\{o} v ' 

the compensated Poisson measure ^ where dt is the Lebesgue measure^. For B G S^IR^-JO}^ where (I B, one 

can show N^it^B) is a martingale and E N u (t,B) = 0, (which stems from v(B) < oo). See Sato [12] . Bertoin [3], 
or Applebaum [5] for further information. 

We are ready to stochastically perturb the systems above. Define W(t) as a standard Brownian motion and a 2 
as the variance of the contacts. Furthermore, take N(dt, dy) as a Poisson measure independent of W(t) and v(-) as 
the intensity measure. We assume that v is a Levy measure and u(M) < oo. Lastly, take h{y) as the affect of random 
jumps in the population, where — 1 < h(y) < 1 for every i/gR, and h(y) is continuously differentiable. For the SIS 
model, Equation (1) becomes the stochastic differential equation 

dS(t) = (- pS(t)I(t) - fiS(t) + m + AJ(t) ) dt - aS{t)I{f)dW(t) + f h(y)S(t-)I(t-)N(dt, dy), 

(3) 

dl{t)= ( j3S(t)I{t) - (A + ii)I(ty\dt + aS(t)I(t)dW{t) - [ h(y)S(t-)I(t-)N(dt,dy). 



For the SIRS dynamic, since the anomaly term captures the effects of quarantine, we will take a function similar 
to h{y), call it j(y), with the same assumptions, except that j(y) is assumed to be nonnegative. Equation (2) is now 
a stochastic differential equation of the form 

dS{t) = (- pS{t)I(t) + SR(t))dt - aS{t)I{t)dW(t), 

dl(t) = (pS(t)I(t) - Mit^dt + aS(t)I(t)dW(t) - Jj(y)I{t-)N(dt,dy), (4) 
dR(t) = (\I(t) - SR{t)^dt+ / j(y)I(t-)N(dt,dy). 



In this paper, we will show that Equation (3) and Equation (4) are well-defined, and we give conditions for 
stability of the disease free equilibrium for both dynamics. 

Given that the initial condition is in the interior of the simplex, we will show that for all finite time, Equations (3) 
and (4) are almost surely in the simplex. We will only show this property for the stochastic SIS, since the stochastic 

SIRS model may be shown in a similar manner. Define K(t) = (S(t),I(t)), A 2 = {y e M 2 : xi,X2 > and x± + x 2 — 
l|, and [•, •] as quadratic variation. For simplicity, we will define S c (t) and I c (t) as the continuous part of the process. 



Proposition 1.1. For all finite t, given that x G A2, P x \K(t) G A2J = 1. 

Proof. Consider the mapping on the simplex G(x) = x\ + X2, and define U(t) — G(L(t)). Ito's lemma yields 

r)C 1 r)C 1 f) 2 C 

dU(t) = —(L(t))dS c (t) + 2 — (L(t))dI c (t) + (£(*)) [dS c (t),dS c (t)] 

+ l^(m)[dlc(t),dl c (t)] + \^^(L(t))[dS e (t),dI e {t)] + lJ^(L(t))[dlS),dS c (t)] 

G(L{t) + (h(y)S(t-)I(t-),-h(y)S(t-)I(t-))) - G(L(t)) 



N(dt, dy) 

= ( - )'.S'(/)/(/) - + + AA/S ).// - aS(l)Iil)d]V(l) + [/3S(t)I(t) - (A + n)I(t)\dt + aS(t)I(t)dW(t) 

N(dt,dy) 



(S(t) + h(y)S(t-)I(t-) + I(t) - h(y)S(t-)I(t-)) - (S(t) + /(*)) 
= ( -Js(t) + I(t)) +fj)dt. 
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Thus, if the process is in A 2 , then dU(t) = ( - (i[S(t) + I(t)j + fijdt = ( - n)dt = 0. Therefore, if L(t) £ A 2 
then U{t) = 1. 

Finally, we show that L(t) does not hit or jump over the boundary in finite time. Define \&(x) = \og(x2/ x\), r as 
the first time L(t) leaves the open simplex (i.e., such that I(t) < or S(t) < ), and Z(t) :— fy(L(t)) for t < r. We 
will apply Theorem 2.1 in Meyn and Tweedie [10] to the process Z{t) in order to show this Z{t) does not explode in 
finite time, and thus P x (t = oo) = 1. By Ito's lemma we have 

dZ (t) = + p-(L(t))dIc(t) + l^(R(t))[dS c (t),dS c (t)] 



dxi 
Id 2 * 



2 dx\ 



(L(t))[dl c (t),dl c {t)] 



dx 2 v ' 2 dx\ 

i a 2 * 



2 dx\dx2 



(L(t))[dS c (t),dI c (t)] 



1 5 2 * 

2 dxidx 



(L(t))[dI c (t),dS c (t)] 



*(i(t) + (h(y)S(t-)J{t-), -h(y)S(t-)I(t-))) - *(L(i)) 



AT(<ft, dy) 



-1 



s . ( / ( ^ - pS(t)I(t) - nS(t) + ^ + XI(t)j dt - aS(t)I(t)dW(t)J 



I(t)\ 



pS{t)Iit) - (A + fi)J(t) (ft + <rS(i)I(t)dW(t) 



-S 2 (t)I 2 (t)dt 



2S 2 (t) w w 2/ 2 (<) 
/(i) - h(y)S(t-)I(t 



S 2 (t)I 2 (t)dt 



log 



S(t) + h(y)S(t-)I(t-) 



- log (/(t)/S(i)) 



N(dt,dy) 



(il{t) + ii+ - X^p- ) eft + o-/(t)dW(t) + ( /3S(t) - (A + /i) ) + aS(t)dW(t) 



^ x {Z{t))+n 



S(t) S(t) 
S 2 (t)dt + 



+ a(V7 1 {Z(t))+^ 1 (Z(t)))dW(t) 



1 + %)/(*-). 



/ l_-%)^(Z(t)) \ 

l0 Hir%j^wr (dMy) ' 

where tf-^s) = _L-(l,e«) = (V (y) , tj" 1 (l/)) • 

Now, defining 2$ as the infinitesimal generator for Z(t) and V(x) = 1 + a; 2 , we see that 

*V{x) = (pHr? (x) + /i + - Xe x + 1 (*) - (A + A*) + y *r' GO* + 1 ') ^' 



a a (*r 1 (*)+^ 1 (*) + 



x + log 



/l-My)*- 1 ^) 



1 + %)^^ 



i/(dy) 



2/3^ 1 (ar) + 2 m- 



-2/i 



- 2Ae 2; + 2/31-2" 1 (x) - 2(A + ^) + cr 2 ^ 1 ( x ) 



-^ r (,)V2/i 0g r-^(i )K%))x. 



Noting that x (a;)| < 1, l*^ 1 ^)! ^ !> and lo S ( ~ 



<log 



1 — min 



%)\ 



+ h(y)^ 1 1 (x) / V 1 + min ' l (2/)/ 

exists positive constants K such that ( BV(x) < KV(x). Therefore P x (r = oo) = 1. 



one can see that there 



□ 
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2 Analysis of the Approximated SIS Model 

Since I(t) = 1 — S(t) we are able to just focus on S(i), which we are able to rewrite as 

dS(t) = (- 0S(t)(l - S(t)) - /iS(t) +fi + X(l- S{t))^Jdt - aS(t)(l - S(t))dW(t) 
+ f h(y)S(t-)(l- S(t-))N(dt,dy) 

J A 

:=a(S(t))dt + -y(S(t))dW(t)+ [ h(y)S(t-)(l - S(t-))N(dt,dy). 

Jr 

Define r £ = inf |< > : S(t) > 1 - e|. 



(5) 



Claim 2.1. For all < x < 1, we have E x 



< oo. 



Proof. We will follow the proof of Theorem 4.2 given in Imhof [BJ. Take L as the infinitesimal generator for our 
process S(t) and define f(y) = e 7 — e lx , where 7 > 0. Now fix an arbitrarily small e > 0. By Dynkin's formula we 
have that 











E Xo 


s(r e AT) 


— e 7 - E 

— e rj Xa 


\f(s(tj)dt 











We will now determine an appropriate upper bound for Lf(s(t)y To adjust for the possibility of the process 
jumping out of the interval, we will consider a x £ [0, 1 — e 2 ]. Hence 



Lf{x) = - 7 (-/3x(l-a;)-/ix + /i + A(l-x))e 73: - 7 - - ^ — 



,-yx _ P ~<( x + h iy) x {}- x )) 



v(dy) 



( 7 CT 2 a; 2 (l-a;)\ . 
= jhx- M - A - 1 { - >-) (1 - x)e^ 



=7z _ p -y(x+h(y)x[l-x)) 



v{dy) 



Noticing that x + h(y)x(l — x) is positive and recalling the inequality — e x < —1 — x for x > 0, we find that 



< 



J2 t- r 



v(dy) < 



1 — ^{x + h(y)x( 



Hence 



■ i=2 

< 7(1 -x)e lx 



Lf{x)< 1 \Bx + 
= j\f3x + 



1 jh(y)x(l-x) 
1 h(y)x 



(l - x) e 72 



v(dy). 



i=2 



1-X 



(7a;) 1 1 h(y)x(l - x) 



v{dy) 



e rx 



v(dy) 



1/(1 


t) 


X 


(1- 


«) 


e jx 







X 


(1- 


x) 


e jx 



h{y)v{dy) - - 1 ^ ^ - M - A (l - x)e^ x 



7cr 2 x 2 (l - x) 



(1 - .T)e 7:E - 7 (/i + A) (l - x)e 7a: . 



i/(M) x /* 7<r 2 x 2 (l x) 

Finally, taking 7 large enough so that f3x + ^ — / h(y)v(dy) < for all x G [0, 1 — e 2 ] 



1-x e 72 



yields 



Lf(x) < -^(ji + xyi-x), 



,-fx 
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Thus 



0< E v 



(t £ At) <e~> -e 2 -/((i + \)E y 



T f AT 



and therefore taking T — > oo, the bounded convergence theorem yields -E Xo r e < oo 



□ 



Theorem 2.1. Suppose that / h(y)v(dy) < 0. If f3 < (j, + \ + J R h(y)v(dy) , then for x e A 2} 



P x lim *T(t) = (1,0) =1. 
\ t— >oo y 

Proof. In our proof we will employ the stochastic Lyapunov method. We will focus on 5(f), defined by Equation(5). 
Define g(x) = 1 — x as our Lyapunov function. Taking Lq as the infinitesimal generator for S(t), we see that 



v{dy) 



L g(x) = ( - /3x(l - z) + - x) + A(l - x)) (-1) + j {l - (s + %)a:(l - x))} - (1 - sc) 
= -^-,3x + ^ + A + y li(y)v((i|/)»j(l-x) 



< 



.{-0 + H + X + J h(y)v(dy)\ (l-x). 



Therefore, Theorem 4 in jlj (page 325) tells us that for an e > 0, there exists a neighborhood of (1,0), say U, such 
that 

PJ lim K(t) = (1,0)) >l-e. 

for x e {/ n A 2 . 

Now, take an arbitrary e > and x e A 2 , and define M = j lim^oo K(t) = (1, 0)|. The strong Markov property 
tells us that 



P X (M) = E* 
Since e was arbitrary, the theorem follows. 



^K(r t ) [XAl] 



> 1 - e. 



□ 



Remark 2.1. Theorem 4 in is stated for a jump-diffusion with a compensated Poisson measure. However, since 
we assumed that i/(K) < oo, we may rewrite Equation (1) as 

dS(t)= ( -pS(t)I(t)-nS(t) + n + \I(t)+ f h(y)S(t-)I(t-)p(dy))dt-aS(t)I(t)dW(t) + f h(y)S(t-)I(t-)N(dt,dy), 
V Jr J Jr 

dI(t)=[PS(t)I{t)-{\ + ii)I(t)- f h(y)S(t-)I(t-)v(dy))dt + aS(t)I(t)dW(t)- f h(y)S(t-)I(t-)N(dt, dy), 



and thus we may apply this theorem. The generator remains unchanged. 



Corollary 2.1. Suppose that h(y) nonnegative, and there exists < ip < 1 such that j3 < \i + A + ip f R h(y)v(dy). 
Then 

PJ lim K(t) = (1,0)) = 1. 

V t— S-OO / 
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Proof. Taking a neighborhood U = {x G A2 : xi > <^j} and as above, we have that 



L g(x) <-(- P + v + X + p h(y)v{dy)\ (l - x) 
Theorem 4 in [4] and Remark 2 gives us, for x € U, 



F x lim tf(t) = (l,0) >l-e , 

1 t— >oo 



for some eg . The rest of the proof follows as above 



□ 



Remark 2.2. The Remark 2 J4H holds for continuous processes. However, since the jumps are positive, this only 
helps the convergence of the sample paths to (1,0). Therefore, we are able to use the result. 



The computer simulations below agree with the conclusion of the theorems. Figure 1 shows for both cases that 
the disease free equilibrium is globally asymptotically stable. Figure 2 tell us that the process is recurrent and thus 
simulates an epidemic. 
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(a) The negative integral case. 



(b) The positive jump function case. 



Figure 1: For simulations, the initial condition arc 5 — -6 and Iq = -4. The parameter values for Figure 1(a) are 
/3 = .1, a = .3, A = .3, jU = .2, v(R) — 1, and j(y) = —.01, while the parameter values for Figure 1(b) are f3 = .4, 
a = .3, A = .3, /i = .1, f(K) = .5, and j(y) = .1. 
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(a) The negative integral case. 



(b) The positive jump function case. 



Figure 2: Just as in Figure 1, the initial condition are So = -6 and io = -4. For Figure 2(a), the parameter values 
are /3 = .4, <r = .3, A = .3, ft = .15, v (RJ = 1, and j(y) = — .1, and for Figure 2(b), the parameter values are f3 = .8, 
u = .3, A = .2, ft = .1, i/(K) = 2, and j(y) = .1. 

Remark 2.3. If / h(y)v(dy) > and i/ie function j(y) was able to take both positive and negative values, we were 
Jr 

not able to say anything about this case. A way to calculate the stability of a system of this type is to determine if 
and where S(t) leaves an arbitrary subinterval of [0, 1]. So, for any < x% < x% < 1, define t XiX2 (xq) — inf |iS(t) ^ 

(xi,xi)\S(0) = x X, ir X2 ; Xl {xo) = P(s(T Xl a; 2 (xo)) >x^j, and n Xl - lX2 (x ) = p(s(t x1X2 (x )) <x^j. The papers of 
\1 6] / and JTjj tell us that solving the integro- differential equation 



a(x)u'(x) + ^^u"(x) + 



u [ x + h(y)x(l — x) ] — u(x) 



v{dy) = 0, 



with the initial conditions of u(x) = 1 for x G [^2)1] o,nd u(x) = for x G [0,:ei], will give us TV Xa ;xi( x o) (and 
accordingly for tt Xi . X2 (xq)). Taking xq — > 0, and x% — > 1 will tell us how this system evolves. This is similar to the 
method given in Gihman and Skorohod J^j. 



3 Analysis of the SIRS Model 

To analyze the SIRS model, we will use the stochastic Lyapunov method. For a bit of simplicity, we will define 
<£(t) = (S(t),I(t),R(t)). 



Note 3.1. Looking closely at the process, one can see that the only stationary position is the point e\ — (1,0,0). 
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Theorem 



3.1. If (3 < mmL + f j{y)v{dy) - j 3 ^v{dy) - s\ then 



P x lim C(t)= ei = 1, 

i->oo 



where i€ A 3 . 



Proof. The proof will follow the one given for the SIS model. For £ the infinitesimal generator of the SIRS stochastic 
process, we see that 



£g(x) = ( - Px-lX 2 + Sx 3 "j t^(x) + {fix x x 2 - \x^j^-(x) + (a.t 2 - 8x 3 ^j^-(x) + \ a2x i x l^2 



(x) 



>d 2 g, 



d 2 g 



, x 2 2 2 U i> l \ 2 2 2 ^ » / \ x 2 2 2 " » / \ 



' dx\dx2 

i(x + (0, -j(y)x 2 , j(y)x 2 )^ - g(x) 



d 2 g 



' dx 2 dxi 



v{dy) 



Take k as a positive constant such that (5-/3 — 2k > 0, (which exists by our assumption). Now, define the positive 
function /(a;) = ci(xi — l) + c 2 x\ + c 3 :rj, where ci, c 2 , and C3 are positive, 



< Cl, 



and 



Then 



(5-/3 -2k 
ci(£ +0) +c 3 / R j 2 (y)^(dy) 

A + / K j(j/)K^)-/ R ^Kdy)-f -P 



< c 2 . 



= ( _ (3xix 2 + 5x 3 ^j2c 3 {x x - l) + ^/3.X!.t 2 - Ax 2 ^c 2 2: 2 + (Xx 2 - Sx 3 ^jc 3 x 3 + \^o 2 x\x 2 2 2c Y + \^r 2 x\x 2 2 2c 2 



+ / c 2 (x 2 - j(y)x 2 ) -c 2 x 2 v(dy)+ f c 3 (x 3 + j(y)x 2 \ - c 3 x 2 3 

JR V 7 . JR V ' 



vidy) 



c\o 2 x\x\ — 2c 2 x 2 + c 2 



2j(y)+3 (y) v{dy)x 2 +c 2 a x x x 2 + 2c 2 j3x 1 x 2 + c 3 / j {y)v{dy)x 



[ ] 2 {y) V {dy)x 2 \ 
Jr ) 



+ 



j - 2cxSx 3 + 2ciSxiX 3 + c 3 Xx 2 x 3 + 2c 3 j j{y)v{dy)x 2 x 3 - 2c 3 (5x 3 j + | - 2cifix 2 x\ + 2ci/3x 2 xi|. 



Organizing and simplifying yields 



£f(x) = 2{c l —a? 1 -c 2 [ A + / jivHdy)- 



v{dy) 



-x x - /3x! + c 3 



[ j 2 (y>(dy)\ 

jr ) 

+ 2 jc^ - 1 + x^ + c 3 (\x 2 + J j{y)v(dy)x 2 - 6x 3 ^j jx 3 + 2ci/?j - x x + l^x 2 x 1 
= 2L^x 2 -c 2 (x+ [j(y)v(dy)-[ j2 Mu(dy)-^x 2 -l3x 1 )+c 3 [ j 2 {y)u{dy)\ 
+ 21 cxS^ - x 2 - x 3 ) + c 3 ( Xx 2 + / j{y)v{dy)x 2 - 8x 3 \ \x 3 + 2ci/3 1 x 2 + x 3 \x 2 Xi 



x\ 



x\ 



2< i ci ( y x * + /3xi ) ~ 02 ( A + / rty) v ( d y) 



Jr z z J jr 



(y)v(dy) \x 



+ 2 j ^ - ci (<5 - ^x ia ; 2 ) + c 3 A + J j(y)u(dy) 



x 2 - (ci + c 3 )5a; 3 ^x 3 



Increasing the values of the positive numbers, we see that 

.2 \ / , rf (y) 



£f(x) < 2 



= 2 



< 2 
- 2 



ci( y + /3 I -c 2 ( A+ / j{y)v{dy) 



"<<///> +C3 / 3 2 (yHdy)\4 



A+ / j(yHdy) 



ei(y+j8j-C2(A+ / j{y)v{dy) 



X 2 - (Cl + C 3 )<5x 3 ) J>X 3 

j 2 (y) 



, y-/ 3 ) + c 3 / j 2 (y)K^)^ 



A+ / j(yHdy) 



x 2 — (ci5 + c 3 (5):r 3 ^ jX3 — kc^x — I) 2 + kc 1 (.e — l) 2 



c i(y + / 3 + K )- c 2(A+/ ]{y)v{dy) - 



(d-P- 2k) 



c 3 



A+ / j{y)v{dy) 



c 1 ^—+f3 + Kj-c 2 ^\ + Jj(y)u(dy) 
c\(5 — n) + c 3 <5^a; 2 — kci(x — l) 2 . 



^Jx 2 ~- (ci (8 — k) + c 3 S^J x^ \x s - kci (x - l) 2 
J 2 (y) 



f (dtf) - y - /3 ) + c 3 / f(y)u{dy) [>:r 



Therefore, there exists a positive constant k, such that £g < —kg. Theorem 4 in [3] tells us that for an e > 0, there 
exists a neighborhood of ei, say U, such that 

P x [ lim <£(i) = ei ) > 1 - e. 

V t— >oo y 

for x € {/ n A 2 . 

Define r e as above, and fix an arbitrary e > 0. We will show that 

adjust for the process possibly jumping over the boundary of r e , we will define I — eo = sup <T J S' (t^) V , (where eo = 1 
if the set is empty). Since the process does not hit the boundary in finite time, we have that e > a - s - Note that, 
for the function / above, Zf(x) < a(l — x) 2 , for some a > 0. Now for x e A 3 , such that x\ < 1 — e, and < t < oo, 
Dynkin's formula yields 



< oo, which will assert our theorem. To 



EL 



/(r e At)J =f{x)-E» 
< f(x) - E» 



zfmt))da 



1-S(t)) ds 



< f(x) - ae 2 E x r e A t 



Thus E x r e < oo, and therefore, following the proof of Theorem 1, 



P x [ lim (S(t) =ei = 1. 

1 t— >oo 



□ 



The simulations given below show globably asymptotic stability to the disease free equilibrium and an epidemic. 
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(a) Convergence to the disease free equilibrium. (b) An example of an epidemic. 

Figure 3: For both simulations, the initial condition are Sq — .3, Iq — .6, and Rq = .1. For the left hand side the 
parameter values j3 — .3, a = .1, A = .29, 5 — .4, = 1, and j(y) = .3 were used. While the parameter values of 

P = .8, a = .2, A = .1, S = .1, v(R) = .5, and j(y) = .1 were used. 
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